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Abstract 



The usual approaches to the definition of energy give an ambiguous result for 
the energy of fields in the radiating regime. We show that for a massless scalar 
field in Minkowski space-time the definition may be rendered unambiguous 
by adding the requirement that the energy cannot increase in retarded time. 
We present a similar theorem for the gravitational field, proved elsewhere, 
which establishes that the Trautman-Bondi energy is the unique (up to a 
multiplicative factor) functional, within a natural class, which is monotonic 
in time for all solutions of the vacuum Einstein equations admitting a smooth 
"piece" of conformal null infinity J7". 
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Consider a Lagrangian theory of fields cf) A defined on a manifold M with a Lagrange 
function density 

C = C[<j> A } «9 M A , ...,d^...d^ A ] , (1) 

for some fceN, where <9 M denotes partial differentiation with respect to x^. Suppose further 
that there exists a function t on M such that M can be decomposed as I x E, where 
£ = {t = 0} is a hypersurface in M and the vector d/dt is tangent to the K. factor. The 
proof of the Noether theorem, as presented e.g. in [Section 10.1] of Ref. |JJ, shows that the 
vector density 

E X — —£X X + X^ J2e=o 4> A ,a!...am 

■<Ey(-iR-^( »*'^ m .., ) < 2 > 

has vanishing divergence, -E a ,a = 0, when the fields 0a are sufficiently smooth and satisfy 
the variational equations associated with a sufficiently smooth £ (cf. also Ref. 0). (This 
is in any case easily seen by calculating the divergence of the right-hand-side of eq. @.) 
Here 4> A tai ...a e = d ai . . . d ae 4> A , and X^d^ = d t . In theories in which £ depends only upon 
<p A and its first derivatives, it is customary to define the total energy associated with the 
hypersurface S by the formula 

E(Z) = J^E x dS x , (3) 

with dS\ —— d\ j dx° A . . . dx 3 , where j denotes contraction. By extrapolation one can also use 
(0) to define an "energy" for higher order theories. Now it is well known that the addition 
to £ of a functional of the form 

d x (Y x [<p A , d a <p A , ...,d ai ...d ak ^ A ]) , (4) 

where k is as in (|TJ), does not affect the field equations^ Such a change of the Lagrange 
function will change E(E) by a boundary integral (see, e.g., Ref. Q3] for an explicit formula 
for A^ A ): 

£(£) — £(E) = + / A^ A d5 MA , (5) 

where S"^ = d a jdp j dx° A ... A da; 3 . If is a "sphere at infinity" the integral over <9£ 
has of course to be understood by a limiting process. Unless the boundary conditions at <9£ 
force all such boundary integrals to give a zero contribution, if one wants to define energy 
using this framework one has to have a criterion for choosing a "best" functional, within the 
class of all functionals obtainable in this way. The vanishing of such boundary integrals will 



1 Here we adopt the standard point of view, that the field equations are obtained by requiring the 
action to be stationary with respect to all compactly supported variations (cf. e.g. Ref. for a 
discussion of problems that might arise when this requirement is not enforced). 
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not occur in several cases of interest, including a massless scalar field and general relativity 
in the radiation regime. 

As an example, consider a scalar field <fi in the Minkowski space-time, with £ = {t = 0}. 
Assume that <fi satisfies the rather strong fall-off conditions on E 

d ai ...d ao <P = o(r- 2 ), < j < A: - 1 , (6) 

where k is the integer appearing in (|l|). In this case the boundary integral in (|5|) will vanish 
for all smooth Y^'s, as considered in eq. (^). This shows that eq. (|3p leads to a well-defined 
notion of energy on this space of fields (whatever the Lagrange function £), as long as the 
volume integral there converges. (That will be the case if, e.g., C has no linear terms in <p 
and its derivatives.) 

Consider, next, a scalar field in Minkowski space-time, with £ being a hyperboloid, 
t = + x 2 + y 2 + z 2 . Suppose further that in Minkowski coordinates C = \rj^ v d^ 4>d u (J), 
so that the field equations read 

□0 = 0. (7) 

In that case the boundary conditions (|6|) would be incompatible with the asymptotic behavior 
of those solutions of eq. (|7|) which are obtained by evolving compactly supported data 
on {t = 0} (see eq. (|9]) below). Thus, even for scalar fields in Minkowski space-time, a 
supplementary condition singling out a preferred E x is needed in the radiation regime. 

Now for various field theories on the Minkowski background, including the scalar field, 
one can impose some further conditions on E x which render it unique HH , such as Lorentz 
covariance, and dependence upon the first derivatives of the field only. The extension of that 
analysis to the gravitational field has been carried out in || , and it also leads to a unique E x 
(namely the one obtained from the so-called "Einstein energy-momentum pseudo-tensor"), 
within the class of objects considered. While this is certainly an interesting observation, 
the restrictions imposed in that last paper on the energy-momentum pseudo-tensor of the 
gravitational field are much more restrictive than is desirable. Thus it seems of interest 
to find a more natural criterion, which would encompass both general relativity and field 
theories on a Minkowski background, and which would single out a preferred expression for 
energy in the radiation regime. In this letter we wish to point out that the requirement of 
monotonicity of energy in retarded time allows one to single out an energy expression in a 
unique way, within a natural class of "energies". Let us start with the case of a massless 
scalar field in Minkowski space-time. The variational formalism described above leads one 
to consider functionals of the form 

Hi<p,t]=E(Z t )+ j H^dS a(3 , (8) 

E{Et) = j T%X v dS^ X v d v = d t , 

St 

where H alS is a twice continuously different iable function of 0(x), <9 ai 0(x), . . ., 
d ai . . . d an (p(x), for some n. The indices a refer to Minkowski coordinates, and in Minkowski 
coordinates the H a/3, s depend upon the coordinates through the fields only. Here and 
throughout the S t 's are unit hyperboloids in Minkowski space-time: = {a; = t + 
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^/(x 1 ) 2 + (x 2 ) 2 + (a; 3 ) 2 }, T* v is the standard energy-momentum tensor for the scalar field 

(with the notation determined by e q . ®), / H°*AS aS is understood as a limit a s R 

est 

tends to infinity of integrals on coordinate balls of radius R included in S 4 . 

Before analyzing convergence of functionals ([8]) we need to specify the class of fields of 
interest. Consider solutions of (^) which have smooth compactly supported initial data on 
the hyperplane {x° = 0}, where x° is a standard Minkowski coordinate. Using conformal 
covariance of eq. (|7|) (cf., e.g., arguments in Ref. J7j) it can be shown that there will exist 
smooth functions c(u, 8, <fi) and d(u, 9, <fr) defined on (— oo, oo) x S 2 such that 

(f)(u, r, 6, 0) = 0(r ), (9) 

with u = x° — r. Moreover (||) is preserved under differentiation in the obvious way. (The 
hypothesis that the initial data are compactly supported is not necessary, and is made only 
to avoid unnecessary technical discussions.) In what follows we will only consider solutions 
of (|7|) satisfying @. For our purposes it is important to emphasize that given arbitrary 
functions c on, say [uq — 1, Uq + 1] x S 2 and do on S 2 there exists a solution of the wave 
equation defined on Minkowski space-time such that (j^) holds, with 

c)d 1 

— = --A 2 c, d(u ,6,4>) = d (9,4>) . (10) 



Here A 2 denotes the Laplace operator on S with the standard round metric. (Eq. ( |10|) is 
obtained by inserting the expansion (^) in ([?])). This assertion is easily proved using again 
the conformal covariance of eq. (0). We claim the following: 

Theorem 1 Let H be as described above and suppose that (W converges for all solutions of 
the wave equation satisfying (j^j. If — ^- < ; then for all such 0's 

H af3 dS af3 = , 

as t 

so that the numerical value of H equals the standard canonical energy. 

PROOF: We can Taylor expand H a ^ at = up to second order to obtain 



H af3 = H af3 \ + V 
1 0=0 *-f 

0<\I\<k 



<Pl 

0=0 



f)2 zia.fi | 

o<i^i<^ 0/90jU=o 

where we use the symbol <pj to denote objects of the form d ai . . . d ae (p, with |/| = 
|(ai, . . . , oti)\ = i. By hypothesis H afS L depends only upon the metric and its derivatives, 

so in Minkowski coordinates the coefficients H a ^ 



0=0 ' ^ 



etc. in (O) are constants. 

0=0 
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By well known properties of Taylor expansions and by (^|) we have r al3 = o(r -2 ), so that r al3 
will not contribute to H in the limit r — > oo. In Minkowski coordinate systems (||) can be 
rewritten as 



C U) 



dx a l ■ ■ ■ dx a i^ ~ r '""i ■ ■ - n aj + \l 1 + 0(r 3 ), (12) 

di^T ■ ■ ■ di^T^ = ^r n ai ■ ■ -n aj + 0(r~ 3 ) , (13) 

where c^(t,r,9,(p) = £^ (c(u, 9, </?)) \ u=t _ r , similarly for S- m \ and n M = (1, — y). Here 
L ai ... a . is a linear functional of c and its tt derivatives up to order j — 1. Inserting (|9|) in (|TTD 
and making use of (^2])-(P^) might produce several terms which do not obviously converge, 
but those have to cancel out or integrate out to zero by our hypothesis of convergence of 
(B). It then follows that (|8|) can be rewritten as 

H = E(E t ) 

+ [ h[c,c^,...,c^,d^,...,d^,9,(f)]d 2 Li , (14) 



! S 2 

for some functional h, smooth in all its arguments, with d 2 [i = sin^d^d^. Moreover h is 
linear in d and its ^-derivatives. Eq. (|H]) allows one to eliminate the w-derivatives of d in 



terms of derivatives of c, so that (|14|) can be rewritten as 
H = E(E t ) 

+ f {h[c,c^,...,c^,9,<p} + a[d,9,<p])d 2 LL , (15) 
Js 2 

for some functionals h and a, with a linear in d. The ^-derivative of ( |I~5"D gives 

d JL = ( (- {c m r + p» + ... 

dt J S 2 \ be 



+ <5c« + 5d ' ^ 



Since c( k+1 > is an arbitrary function on S 2 at fixed c, . . . , c^ k ' and do, we can choose it so 
that ^ < unless Sh/Sc^ — 0, k > 1. A suitable redefinition of ft, leads to 

#(*) = S(£ t ) + J s2 (h[c, 9, 0] + a[d, 9, 0]) d 2 /i , 
f = / S2 (H (1) + t)^ (1) -|iA 2 c)^, (16) 

Consider now solutions of the wave equation with c^'(u = u ) = 0. In this case and 
arbitrariness of c(u = uq) imply that dH/dt will be non-positive if and only if A2(5a/5d) = 
, which forces 5a/ 5d to be a constant. We note that for any constant a the integral 

add 2 fx, (17) 

s 2 

is a constant of motion (see also || Sect. 8.2]. However, integrals of the form (|i~7|) cannot 
arise in the class of functionals considered here. Indeed, the identity (|H|) shows that all the 
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terms which would give a non-vanishing contribution to H and which contain derivatives of 
d contain at least one u derivative of d. Then the only possible term which would contain d 
would come from the term 

f dHaP j. ja f dRrt / * J2 

/ fl i 9 dS <*0 = / -rrr {cr + d)d // , 

J S 2 d(f) <f>=o J S 2 d(j) 0=0 

which for generic c diverges when r goes to infinity, unless identically vanishing. We thus 
obtain 5a/5d = 0. The right hand side of eq. ( |l6|) can be made positive by choosing 
c^(u = uq) = ^Sh/Sc, unless Sh/Sc = 0, and our claim follows. □ 
Let us now turn our attention to those gravitational fields which are asymptotically 
Minkowskian in light-like directions. An appropriate mathematical framework here is that 
of space-like hypersurfaces which intersect the future null infinity Sf + in a compact cross- 
section K. For such field configurations it is widely accepted that the "correct" definition of 



energy of a gravitating system is that given by Freud M, Trautman |10|,|TT|, Bondi et al. [12 



and Sachs [I3|, which henceforth will be called the Trautman-Bondi (TB) energy. There 



have been various attempts to exhibit a privileged role of that expression as compared with 
many alternative ones ( |I^-p^,|^,p^-p3[] , to quote a few), but the papers known to us have 
failed, for reasons sometimes closely related to the ones described above, to give a completely 
unambiguous prescription about how to define energy at $. (We make some more comments 
about that in M, cf. also p3fl .) In a way rather similar to that for the case of the scalar 
field described above, in ||] we show that the TB energy is, up to a multiplicative constant 
«6l, the only functional of the gravitational field, in a certain natural class of junctionals, 
which is monotonic in time for all vacuum field configurations which admit (a piece of) a 
smooth null infinity £[ + . More precisely, in || we show the following: 

Theorem 2 Let H be a functional of the form 

H[g, u}= H a(S (g^ u , g mtr , ^ j(J1 ... CTfe ) dS Q(3 , (18) 

JS 2 (u) 

where the H a/3 are twice differentiable functions of their arguments, and the integral over 
S 2 {u) is understood as a limit as p goes to infinity of integrals over the spheres t = u + p,r = 
p. Suppose that H is monotonic in u for all vacuum metrics g^ iu for which H is finite, 
provided that g^ u satisfies 



-2\ 



d n ■ ■ ■ d ai {9»» - "-^y^ ~ n -^^) = o(r" 2 ), (19) 

with 1 < i < k, for some C k functions h a ^ u (u, 6, <p), a = 1,2. If H is invariant under 
passive BMS super-translations, then the numerical value of H equals (up to a proportionality 
constant) the Trautman-Bondi mass. 

Some comments are in order. First, the volume integral E(H t ) which was present in 
(^) does not occur in (0), because the Trautman-Bondi mass is itself a boundary integral. 
Next, Theorem imposes the further requirement of passive BMS invariance, which did not 
occur in the scalar field case. This requirement arises as follows: recall that the coordinate 
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systems in which the metric satisfies (|T9D are, roughly speaking, defined only modulo BMS 
transformations. Then the requirement of passive BMS invariance is the rather reasonable 
requirement that the concept of energy be independent of the coordinate system chosen to 
measure this energy. We note that we believe that the requirement of monotonicity forces 
the energy to be invariant under (passive) super-translations, but we have not succeeded in 
proving this so far. 

The proof of Theorem ^] is similar to the proof for the scalar field presented here, but 
technically rather more involved. A key ingredient of the proof is the Friedrich-Kannar 
P5|-P7| construction of space-times "having a piece of JT' . 

It is natural to ask why the Newman-Penrose constants of motion |28| , or the logarithmic 
constants of motion of , do not occur in our results of [|J . These quantities are excluded by 
the hypothesis that the boundary integrand H a @ which appears in the integrals we consider 
depends on the coordinates only through the fields. The Newman-Penrose constants could 
be obtained as integrals of the form (cf. e.g. [^U|) if explicit r 2 factors were allowed in 
H af3 . Similarly logarithmic constants could occur as integrals of the form (|2]) if explicit 
1/ In r or r +l In - - 7 r factors were allowed there. 

Let us finally mention that one can set up a Hamiltonian framework in a phase space 
which consists of Cauchy data on hyperboloids together with values of the fields on appro- 
priate parts of Scri to describe the dynamics in the radiation regime [31]. Unsurprisingly, the 
Hamiltonians one obtains in such a formalism are again not unique, but the non-uniqueness 
can be controlled in a very precise way. The Trautman-Bondi mass turns out to be a Hamil- 
tonian, and an appropriate version of the uniqueness Theorem || can be used to single out the 
TB mass amongst the family of all possible Hamiltonians. In the Hamiltonian framework 
the freedom of multiplying the functional by a constant disappears. 
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